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I. INTRODUCTION 


The electroweak-scale right-handed neutrino (EW-scale uji) model was proposed 
by one of us (PQH) [T] with the following main motivations in mind: 1) To provide 
a model for the see-saw mechanism which can be realized at the electroweak scale 
instead of a typical grand unihcation theory (GUT) scale; 2) To be able to test the 
seesaw mechanism through the discovery of right-handed neutrinos whose Majorana 
masses are now bounded by the electroweak scale Aew ~ 246 GeV; 3) To be able 
to probe at high energies (e.g. at the Large Hadron Gollider (LHG)) lepton-number 
violating processes such as like-sign dilepton events coming from the Majorana nature 
of the right-handed neutrinos. The electroweak-scale right-handed neutrinos belong 
to doublets of the Standard Model (SM) SU{2) whose partners are right-handed 
“heavy” mirror charged leptons. The requirement of the absence of anomaly dictates 
the addition of right-handed doublets of mirror quarks to the particle spectrum. 
Furthermore, left-handed S't/(2)-singlet mirror quarks and mirror charged leptons 
will be the counterparts of their SM right-handed SU (2)-singlet quarks and charged 
leptons. 

The EW-scale pr model entails extra SU{2) chiral doublets (the mirror fermions) 
which have many consequences. These mirror fermions enter loop corrections to 
various quantities and processes such as the electroweak precision parameters, rare 
processes, etc. 

The hrst type of effects that needs to be examined is the contributions of these 
extra chiral doublets to the electroweak precision parameters. These calculations 
have been performed in [2] and it was found that there is a large parameter space 
where the EW-scale ur model satishes the EW precision constraints. In a nutshell, 
the contributions from the mirror fermions are partially cancelled by those of the 
scalar sector, in particular the SU (2) triplet scalar. 
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The next place where mirror fermions enter through loop corrections is rare pro¬ 
cesses such as /r —)■ 6 7 and r —>■ /i 7 . In |3], such processes have been discussed in a 
generic fashion, with an emphasis on the possible correlation between the observabil¬ 
ity of the aforementioned rare processes and the decay lengths of the mirror charged 
leptons, both of which are of phenomenological interests. In this article, we will 
present an update of the process /x —ey taking into account recent developments 
of the model, including experimental inputs from the recently-discovered 125 GeV 
SM-like scalar BIE]- They are summarized below. 

The scalar sector of the original model pQ contains one SM-like Higgs doublet and 
two Higgs triplets, one with Yj2 = \ containing doubly-charged scalars and one with 
Y/2 = {]. (The rationale for this sector will be explained in the summary section.) 
The discovery of the 125 GeV SM-like scalar has opened up a whole new chapter 
on any model beyond the SM, in particular those models which have more than one 
Higgs doublet. In light of this discovery, a close examination of the scalar sector of 
the EW-scale model [B] revealed that its original Higgs content is insufficient to 
accommodate the 125 GeV SM-like scalar. It turns out that a simple introduction of 
an extra Higgs doublet, now totaling two: one of which couples to the SM fermions 
and the other one to the mirror quarks and charged leptons. This yields two 125- 
GeV candidates with one being SM-like (dubbed Dr. Jekyll) and the other being very 
different [Mr. Hyde), both of which giving comparable signal strengths in agreement 
with ATLAS and GMS data. 

Most importantly for the present manuscript is the recent work [7j concerning 
neutrino and SM charged lepton masses and mixings. The fact that the SM lep¬ 
ton mixing matrix Vpmns (the Pontecorvo-Maki-Nakagawa-Sakata mixing matrix) 
is so different from the quark counterpart, Vckm (the Gabibbo-Kobayashi-Maskawa 
mixing matrix), has given rise to many models, many of which invoke the pres¬ 
ence of some kind of discrete symmetry. Among these different proposals for the 
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discrete symmetry is the popular A 4 symmetry which has been used to reproduce 
the tribimaximal form of t/pMNS- This symmetry is usually applied to the charged 
lepton sector [S] and involves four or more Higgs doublets. (Such a large number 
of Higgs doublets might be hard to accommodate the 125 GeV SM-like scalar with 
the desired observed properties.) The new twist of [7j is to exhibit the A 4 symme¬ 
try in the neutrino Dirac mass sector and the scalar sector involved is composed of 
SU{2) X f/(l)y-singlet scalars which are not constrained by LHC data. These singlet 
scalars are composed of a singlet and a triplet of A 4 . This model reproduces the 
desired PMNS matrix and makes predictions on the charged lepton mass matrix in 
the form of The singlet scalars play a crucial role in the process /r —)■ ey 

in the EW-scale i/r model as shown in j3] and updated below in light the aforemen¬ 
tioned developments. The results presented in this paper contain a deep correlation 
between the branching ratio —)■ ey) and the neutrino sector in the form of the 
PMNS matrix for both normal and inverted hierarchies, as well as the form of the 
mirror lepton mixing matrix. It will be shown that the exclusion zones in the plots 
of the branching ratio of B(/u ej) versus the Yukawa coupling strengths to the 
singlets depend a bit on how strong the y44-triplet scalars couple to the leptons. 

This paper is organized as follows. First, in section H, we summarize the essence 
of the EW-scale model (original [1] and extended 0 ). Next, in section HI, we 
briefly review constraints from electroweak precision measurements for the original 
model and from Higgs physics for the extended model. In section IV, we briefly 
review the results of neutrino and charged lepton masses and mixings [7]. We then 
proceed with the actual calculations of the process —)■ Ijj, the anomalous magnetic 
dipole moment Aa;^ and the electric dipole moment di. for the lepton U in section 
V. Detailed numerical analysis will be presented in section VI. Implications of our 
results concerning the possible detection of mirror leptons at the LHC and the ILC 
are discussed in section VH. We hnally summarize and conclude in section VIII. A 
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few useful formulas are collected in an Appendix. 


II. REVIEW OF THE EW-SCALE ur MODEL 

For the sake of clarity, we review in this section the original EW-scale vr model 
[T] and its extended version p. 

• Gauge group: 


SU{?,)c X SU{2) X U{1)y 


( 1 ) 


There are many differences between the EW-scale pr model and the popular 
Left-Right symmetric model [S|. The hrst difference lies in the gauge structure 
of the two models: SU{3)c x SU{2) x U{1)y for the EW-scale pr model and 
SU{3)c X SU{2)l X SU{2)r x U{1)b-l for the Left-Right model. 

• Lepton and quark SU{2) doublets (the superscript M refer to mirror 
fermions): 




; Mirror: 



• Lepton and quark SU{2) singlets: 


SM: cr] Ur, dR, Mirror: ef; nf, df. 


• Doublet Higgs fields: 
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As explained in [T], a Higgs doublet is needed to give masses to all charged 
fermions. The analysis of the properties of the 125-GeV SM-like scalar necessi¬ 
tates the introduction of one extra Higgs doublet as explained in [B]. Each Higgs 
doublet couples to a different sector: $2 = (0^,02) to the SM fermions and 
^ 2 M = (02M5 02m) to the mirror fermions. They develop the following vacuum- 
expectation-values (VEV): ($ 2 ) = (0,n2/\/2)^ and {^ 2 m) = (0,n2M/\/2)^- 

• Triplet Higgs fields: 

The SU (2)-triplet Higgs helds form the cornerstone of the EW-scale vr model. 
As shown in [T|, the VEV of the V/2 = 1 triplet gives an electroweak-scale 
Majorana mass to the right-handed neutrinos and the V/2 = 0 triplet is needed 
to preserve the custodial symmetry so that the p parameter equals unity at tree 
level. This is summarized below. Here we just write down the triplet Higgs 
helds and their VEVs. 

- X (>"/2 = 1) = ^ T.x = i ^ I with (x°) = vm- 

\ ^ -75^^/ 

- ^ (V/2 = 0) = ■C”) (ill order to restore Custodial Symmetry) with 

{e°) = VM. 

- VEVs: 

-t- Svj, = (246 GeV)T 

• Singlet Higgs fields: 

The original model which is basically concerned with the energy scales which 
enter the seesaw mechanism contains only one SU{2) x f/(l)y-singlet Higgs 
held 05 whose VEV {(ps) = vs gives the Dirac mass to the neutrinos (to be 
summarized below). It was almost a “toy model” in that it did not discuss 
lepton mixings and, in particular, the PMNS matrix Vpmns- This problem 
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has been recently investigated by [7] within the framework of an A 4 symmetry 
which is applied to the neutrino sector of the EW-scale model. The upshot 
of this study was the introduction of an A 4 singlet 0 o 5 and an y44-triplet 
(i = 1,2,3). They obtain the following VEVs vq and Vi respectively. We will 
summarize below the main points concerning this singlet scalar sector in the 
construction of tlpMNS and its implication to rare progresses such as /r —>■ ey. 

• Dirac neutrino mass 

For simplicity, we will denote the right-handed neutrino helds by i/ji from 
hereon. 

The original model contains one singlet scalar whose VEV provides a Dirac 
mass for the neutrinos. A generic Yukawa coupling is of the form 


= —Qsi II 4>s Ir + H.c. 

= -gsii^L + hi e^) 05 -h H.c. 


( 2 ) 


With (05) = Vs, one obtains the Dirac mass = gsi vs- 

• Majorana neutrino mass 

This is the main point of [1]. The electroweak-scale Majorana mass for the 
right-handed neutrinos is obtained from the following coupling 



(3) 


With (x°) = vm, the Majorana mass is given by Mr = Qm^m- 
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III. REVIEW OF RESULTS OF THE EW-SCALE ur MODEL AS DIS¬ 
CUSSED IN [ 2 \ AND |6| 


In this review section, we will discnss two sets of results for the EW-scale ur model 
obtained in [2] (the electroweak precision constraints) and [6] (constraints from the 
125-GeV SM-like scalar). 


A. Electroweak precision constraints on the EW-scale ur model |2j 


The presence of mirror quark and lepton S'f/(2)-doublets can, by themselves, 
seriously affect the constraints coming from electroweak precision data. As noticed 
in |2], the positive contribution to the S-parameter coming from the extra right- 
handed mirror quark and lepton doublets could be partially cancelled by the negative 
contribution coming from the triplet Higgs helds. Ref. |2] has carried out a detailed 
analysis of the electroweak precision parameters S and T and found that there is a 
large parameter space in the model which satishes the present constraints and that 
there is no fine tuning due to the large size of the allowed parameter space. It is 
beyond the scope of the paper to show more details here but a representative plot 
would be helpful. Fig. 1 shows the contribution of the scalar sector versus that of 
the mirror fermions to the S-parameter within la and 2cr. In this plot, [2] took for 
illustrative purpose 3500 data points that fall inside the 2a blue region with about 
100 data points falling inside the Icr red region. More details can be found in [2]. 


• 1o constraint 
: 2 o constraint 



FIG. 1. Constrained Ss versus Smf- 


B. Review of the scalar sector of the EW-scale vr model in light of the 
discovery of the 125-GeV SM-like scalar [6| 


In light of the discovery of the 125-GeV SM-like scalar, it is imperative that any 
model beyond the SM (BSM) shows a scalar spectrum that contains at least one 
Higgs held with the desired properties as required by experiment. The present data 
from CMS and ATLAS only show signal strengths that are compatible with the SM 
Higgs boson. The dehnition of a signal strength /x is as follows 


/i(iL-decay) 


cT(iL-decay) 
crsM(^f-decay) ’ 


( 4 ) 


with 


cT(iL-decay) = (T(iL-production) x H(iL-decay). 


( 5 ) 


To really distinguish the SM Higgs held from its impostor, it is necessary to 
measure the partial decay widths and the various branching ratios. In the present 
absence of such quantities, the best one can do is to present cases which are consistent 
with the experimental signal strengths. This is what was carried out in [B]. 
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The minimization of the potential containing the scalars shown above breaks its 
global symmetry SU{2)l x SU{2)r down to a custodial symmetry SU{2)£) which 
guarantees at tree level p = M^/M‘^cos^9w = 1 |S]- The physical scalars can be 
grouped, based on their transformation properties under SU ( 2 )o as follows: 

five-plet (quintet) —)■ Hf, 

triplet —)■ 
triplet ^ 

three singlets —)■ ■ ( 6 ) 

The three custodial singlets are the CP-even states, one combination of which can 
be the 125-GeV scalar. In terms of the original helds, one has — 02 M 

and H^' = j. These states mix through a mass matrix obtained 

from the potential and the mass eigenstates are denoted by if, H' and if", with the 
convention that the lightest of the three is denoted by if, the next heavier one by 
H' and the heaviest state by H". 

To compute the signal strengths p, Ref. | 6 ] considers H —)■ ZZ, W~^W~, 77 , bb 
and rf. In addition, the cross section of gg ^ H related to H ^ gg was also 
calculated. A scan over the parameter space of the model yielded two interesting 
scenarios for the 125-GeV scalar: 1) Dr. JekylFs scenario in which H ~ meaning 
that the SM-like component = cf)^ is dominant] 2) Mr. Hydds scenario in which 
H ~ H^' meaning that the SM-like component if® = (j)^^ is subdominant. Both 
scenarios give signal strengths compatible with experimental data as shown below in 
Fig. (2). 

As we can see from Fig. (2), both SM-like scenario {Dr. Jekyll) and the more 
interesting scenario which is very unlike the SM {Mr. Hyde) agree with experiment. 
As stressed in [B] , present data cannot tell whether or not the 125-GeV scalar is truly 
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H^bB 

CMS; = 0.93 ±0.49 


CMS:n = 0.91±0.27 

H ^ YY 

CMS;m= 1.1310.24 


H ^ W"W 

CMS;|i = 0.8310.21 


H^ZZ 

CMS: n= 1.0010.29 



0.5 


1.5 2 2.5 

Best fit o / Oqm 


CMS preliminary 
'mu = 125.7 GeV 


EWVp "Dr. Jekyll" Ex. 1 
' mj:j= 125.7 GeV 

EWVr "Dr. Jekyll" Ex. 2 
' m- = 125.7 GeV 


EWVp "Mr. Hyde" Ex. 1 
m~ = 125.8 GeV 

EWVp "Mr. Hyde" Ex. 2 
m- = 125.2 GeV 


EWVp, "Mr. Hyde" Ex. 3 
= 125.6 GeV 


FIG. 2. Predictions of signal strength /i(i? —)■ bb, rf, 77 , W^W~, ZZ) in the EW-scale 
model for examples 1 and 2 in Dr. Jekyll and example 1, 2 and 3 in Mr. Hyde scenarios 
as discussed in [ 6 ], in comparison with corresponding best fit values by CMS [TOUT^ . 

SM-like or even if it has a dominant SM-like component. It has also been stressed 
in [ 6 ] that it is essential to measure the partial decay widths of the 125-GeV scalar 
to truly reveal its nature. Last but not least, in both scenarios, 
subdominant but is essential to obtain the agreement with the data as shown in [B]. 

As discussed in detail in [ 6 ] , for proper vacuum alignment, the potential contains 
a term proportional to A 5 (Eq. (32) of | 6 ]) and it is this term that prevents the 
appearance of Nambu-Goldstone (NG) bosons in the model. The would-be NG 
bosons acquire a mass proportional to A 5 . 

An analysis of GP-odd scalar states H^, the heavy GP-even states H', H" 

was presented in |H]. The phenomenology of charged scalars including the doubly- 
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charged ones was also discussed in [H]. 

The phenomenology of mirror quarks and leptons was briefly discussed in [2] and 
a detailed analysis of mirror quarks will be presented in [Tn|. It suffices to mention 
here that mirror fermions decay into SM fermions through the process —)■ q(j)s, 

—)■ l(j)s with (ps “appearing” as missing energy in the detector. Furthermore, the 
decay of mirror fermions into SM ones can happen outside the beam pipe and inside 
the silicon vertex detector. Searches for non-SM fermions do not apply in this case. 
It is beyond the scope of the paper to discuss these details here. 


IV. REVIEW OF NEUTRINO AND CHARGED LEPTON MASSES AND 
MIXINGS IN THE EW-SCALE i^r MODEL 


Since the ideas and notations coming out of this review will be important for 
the calculation of the rate of /i —>■ ey, we will present a little more details than the 
previous section. 

In 1^, a model of the Dirac part of neutrino masses was constructed using the 
widely popular A 4 symmetry. Unlike previous works on that symmetry where there 
was a need to introduce several (more than two and typically four or five) Higgs 
doublets (see the review by m and where it might be very problematic with the 
discovery of the 125-GeV SM-like scalar, the main motivation of is to first obtain 
the Cabibbo-Wolfenstein matrix [16] 


U 1 U 


1 




V 


(7) 


1 u 

1 U J 

which is a prototype of the PMNS matrix with “large” mixing parameters and which, 
upon a slight modification, could reproduce the “experimental” Upmns being defined 


as 
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UpMNS = UlU[. ( 8 ) 

Under and transform as 3, where e and u are generic 

notations for the charged and neutral leptons. Using the A 4 multiplication rule 
3 X 3 = 1(11 + 22 + 33) + i'(ll + u^22 + u33) + 1"(11 + a;22 + u^33) + 3(23,31,12) + 
3(32,13, 21) with u = it was argued in that the appropriate set of singlet 

scalars is composed of an ^4 singlet (pos and an y 44 -triplet {(pis} {i = 1,2,3). To 
reflect the two different ways that the 744 -triplet can couple to the leptons, [7j wrote 
down the Lagrangian 

^ { 9 os<Pos + 9 is^s + 92 s^s) Ir’^ + H.c., (9) 

where and are gauge eigenstates which are related to the mass eigenstates by 

ll = Uih . iV = U^rI^. ( 10 ) 


Using the aforementioned multiplication rule, one obtains the following matrix 


9os4>OS 9is4’3S 92s4>2S ^ 

= 92S<p3S 9os<Pos 9is<Pis 

y 9lS(p2S 92S(plS 90S<p0S 

As shown in [7j, reality of neutrino Dirac masses implies that 

92 s = 9is ■ 


( 11 ) 


( 12 ) 


Furthermore, it was shown that, with vq = {(pos) and Vi = {(pis) = v, the neutrino 
mass matrix 


mF = 


^ 903^0 9isVo 92SV2 ^ 
92S'V3 9oS'Vo 913^1 
9iS'V 2 923^1 9o3'Vo 


(13) 
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can be diagonalized, i.e. UlM^Uy, by the matrix 


/i 1 




1 cn 
1 a; 


( 14 ) 


Notice that Ui, = Eqs. (14) and (11) will form a basis for onr snbseqnent 

discnssion. 

For the pnrpose of the snbseqnent sections, we rewrite Eq. 0 as follows 

Cs = -‘l + H.c. (15) 

= -h + Hx., (16) 


where 


and 


= UlM^U, , 


(17) 


~ ^IUr ■ (18) 

The above constrnction can be straightforwardly generalized for the right-handed 
leptons and left-handed mirror leptons. Hence the total £5 becomes 

£s = —II h^PMNS ^4) ~ ^PMNS ^4> ^PMNS^L + (19) 


where M'^ = UlM’^Uu and is the same as given by Eq. (11) with gos —^ g'osy 
gis g'ls S'ss g2s- Reality of the eigenvalnes of also implies 5^25 = 5 ' 15 - 
In analogous to t/pMNS and we have dehned the following mixing matrices 


for the second term of Eq. (19) 


and 


^PMNS ~ ^l^Ri 


tt/M _ TjUjl^ 
'-^PMNS ~ L 1 


( 20 ) 


( 21 ) 
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where and f/|f^ are the unitary matrices relating the gauge eigenstates and the 
mass eigenstates 


I 


0 _ 
R — 


KIr 


iMfl = 


( 22 ) 



FIG. 3. One-loop induced Feynman diagram for h —)• in EW-scale vr model. 


V. THE CALCULATION 


The one-loop irreducible diagram for U —)■ is shown in Fig. ([^. Other two 

diagrams not shown are reducible associated with the one-loop dressing for the ex¬ 
ternal fermion lines. They are crucial for the cancellation of ultraviolet divergences 
and gauge invariance in our calculation. The relevant Yukawa couplings between 
the leptons, mirror leptons and the singlet and triplet scalars can be deduced by 


recasting the Lagrangian Cs in Eq. (19) into the following component form 

3 3 

= ('‘k + UiUSl'L) ‘t’ts + H.C. 

k=0 i,m=l 


(23) 
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where 


— (^PMNS ■ ■ ^PMNs) ) (24) 

\ / im 

3 

^ (^PMNs) . . (^™Ns)nm ’ (^^^ 

T ' '^3 

J,n=l 

and 

= (^PMNS ' M ^ ■ t/pMNs) ) (26) 

= X] (^PMNs) .^. ^jn (^PMNs)nm ' (^'^) 

j,n=l 

The matrix elements for the four matrices M^{k = 0,1, 2, 3) are listed in Table I. 
Mj^ can be obtained from Mj^ listed in Table I with the following substitutions 
90S 9os and gis ^ g'ls- 


A. The process h —)■ Ij'y {i ^ j) in EW-scale vr Model 

Lorentz and gauge invariance dictate the form of the amplitude for the process 

C(P) ^ ^J{P')+^{9) to be 

M (/■ ^ IJ 7 ) = €*^{q)uj{p) {ia^^'^qu [C"IPl + C'^Pr] } Ui{p) , (28) 
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where Pl,r = (1 =F 75)/2. The coefficients C^l^ can be extracted from the one-loop 
diagram (Fig. 


C'l = 


EE 


Ihvr^ ' I 

k=0 m=l I 


m,ufj (US„^y + («“)' 


m 


4‘kS 


m 


[M 




m 


<f>kS 


c-> = + 


EE 


1 


m 


(29) 


[M 


Ihvr^ ^' I m,M '- 

k=0 m=l I 


7“)' + («“)' 


m 


</'fcS 


m 


IM 




m 


</>fcS 


rriiM 

‘'m 


m 


(30) 


/M 

‘'m 


Here we have assumed the mirror lepton masses are much larger than the external 
fermion masses rriiM rriij and set rriij —)■ 0 in the loop functions X(r) and J7’(r), 
which are simply given by 


I(r) = 


[—6r^ logr -|- r(2r^ -|- 3r — 6) -|- l] , 


12(1 -r)4 

(^) ^ n/i ^ [~2r^ log r -I- r(3r — 4) -I- l] 


(31) 

(32) 


2(1 -r)3 

In our numerical work for /r —)■ ey presented in section VI, we will consider the 
mirror lepton masses of the order a few hundred GeV and the singlet and triplet 
scalar masses of the order 10 MeV, thus the ratio r = ~ 10“^ is very 


tiny. For all practical purposes, one can replace Eqs. (31) and (32) by the limits 


lim,._,.oX(r) = 1/12 and Xmir^oJir) = 1/2 respectively. Formulas of X and J for 
the general case of rriij ^ 0 are given in the Appendix. 

The partial width for /j —)■ Ij'y is given by 

3 


F {k —)• Ij'y) — 


m. 


IOtt 


-m. 


l-Jr (K'P + ICgP) . 


(33) 
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B. Magnetic Dipole Moment 


The magnetic dipole moment anomaly for lepton /j can be easily extracted from 
the above calculation with the following result 

2m,, (Cl + Cf 


Aai. = 


= d 


2 

fc=0 m=l 


1 


IM \ f^lM 


3 3 


jLk (TjRk\*\ n- I ^<t>ks 


Re J 


k=0 m=l 


rriiM 


m 




(34) 


C. Electric Dipole Moment 


The electric dipole moment operator for a fermion / is usually dehned as 


C 


EDM — , 


(35) 


where is the electromagnetic held strength and the coefficient df the electric 
dipole moment. The electric dipole moment for lepton /j can also be easily extracted 
from the above calculation with the result 




fc=0 m=l 


Lk 


J 


m 


4>kS 


m 




(36) 


VI. NUMERICAL ANALYSIS 


The branching ratio B(/r —)■ ey) is given by 


B(/i ey) = ■ T {k -)■ Ij'j) 


(37) 
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where is the lifetime of the muon im 

= (2.1969811 ± 0.0000022) x 10“® s . (38) 

In our numerical analysis, we will adopt the following approach: 

• For the masses of the singlet scalars (j)ks, we take 

= Ms : 2Ms : 3Ms : AMs 

with a hxed common mass Ms = 10 MeV. As long as rn(j,^g results 

will not be affected much by the exact mass relations among these singlet 

scalars. 

• For the masses of the mirror lepton we take 

~ Afmirror “1“ 

with 5i = 0, ^2 = 10 GeV, ^3 = 20 GeV and vary the common mass M^i rmr 
from 100 GeV to 800 GeV. 

• We assume all the Yukawa couplings goS) 5 'is, 5 ' 25 , g'os^ 9iS) ^e all 

reaQ As mentioned before, g 2 s = {qis)* and g^s = {g'ls)* due to the reality of 
the mass eigenvalues of the Dirac neutrino masses. For simplicity, we also take 
90 S = g'^si 9^s = 9 'is study the following 6 cases: 

1- 90S 7^ 0, gis = 0. The A 4 triplet terms are switched off. 

2. gis = 10“^ X g^s- The A 4 triplet couplings are merely one percent of the 

singlet ones. 

^ In this study, we do not analyze the possibility of electric dipole moments for the charged leptons 
in which complex Yukawa couplings must be assumed. 
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3. gis = 10 ^ X qqs- The triplet couplings are 10 percent of the singlet 
ones. 

4. gis = 0.5 X g^s- The triplet couplings are one half of the singlet ones. 

5. gis = gos- Both ^4 singlet and triplet terms have the same weight. 

6 - 90 S = 0) 9is 7 ^ 0- The ^4 singlet terms are switched off. 

• For the three unknown mixing matrices h^pMNS ^pmnsj 

consider two scenarios: 


— Scenario 1 


TjM _ Tj! _ tt/M _ 

'^PMNS ~ '^PMNS ~ '^PMNS ~ 


'^CW 


— Scenario 2 

= ^PMNS = ^PMNS = ^PMNS 


Recall that the standard parameterization of the PMNS matrix is given 
by 


t4>MNS 


C12C13 ■S12C13 §136 

—512023 — 0125235136**^ O12O23 — 512 5235136**^ 523O13 

512523 — 6126235136**^ —012523 “ 5126235136**^ O23O13 y 


where sij = sin 6 *jj, Ojj = cos% and P = Diag(l, 6 *“*^*/^, 6 *"*^*/^) is the 
Majorana phase matrix. We will ignore the Majorana phases in this 
analysis. 


In Table II we list the Icr range of the mixing parameters as given by the 
recent analysis of global three-neutrino oscillation data in [T^ [T^ . With 
the central values for the mixing parameters given in Table II as inputs, 
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we obtain two possible solutions of the PMNS matrix: 

/ 

NH 




PMNS 


V 


0.8221 0.5484 -0.0518 + 0.1439* 

-0.3879 + 0.07915* 0.6432 + 0.0528* 0.6533 

0.3992 + 0.08984i -0.5283 + 0.05993* 0.7415 


A 


for normal hierarchy, and 

( 


U, 


IH 

PMNS 


V 


0.8218 0.5483 -0.08708 + 0.1281* 

-0.3608 + 0.0719* 0.6467 + 0.04796* 0.6664 

0.4278 + 0.07869* -0.5254 + 0.0525* 0.7293 


A 


/ 


for inverted hierarchy. For each scenario, we consider these two possible 
solutions for the Fpmns- Due to the small differences between these two 
solutions, we expect our results are not too sensitive to the neutrino mass 
hierarchies. 

Limits on i?(/* —)■ ey) from MEG experiment [20] and its projected sensitivity 

m- 


ey) < 5.7 x 10"^^ (90 C.L.)[MEG, 2013] , (39) 

S(p —)■ ey) ~ 4 X 10“^^ [Projected Sensitivity] . (40) 

• Aa^ from E821 experiment |22j : 


Aa^ = a;fP - a™ = 288 (63) (49) x 10 


-11 


(41) 


Since the dominant contributions to the loop amplitude arise from the mass in¬ 
sertion of the internal mirror lepton line in Fig. (|^, only the last terms in Eqs. (29), 


(30) and (34) are signihcant numerically. As long as the current MEG 
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TABLE II. Mixing parameters from global three-neutrino oscillation data taken from [181 

m- 


Mixing Parameters 

Normal Hierarchy 

Inverted Hierarchy 

sin^ 012 

0.308 ±0.017 

0.308 ±0.017 

sin^ 023 

n 4Q7+0.033 
u.nnr_o 023 

U.400_o 031 

sin^ 013 

0 n244+0-0020 

U.UZOn_o 0019 

u.uzn_o 0022 

bj-K 

X Qq+0.38 
l.OU-o 27 

1 +0-29 

-‘-*'^-‘■-0.33 

5m? = — m\ 

(7.54l°|6) X lO-^eV^ 

(7.54l°J6) X lO-^eV^ 

Am^ = — (mf -|- m 2)/2 

(2.43 ± 0.06) X lO-^eV^ 

(2.38 ± 0.06) X lO-^eV^ 


limit (Eq. 


on the branching ratio i?(/i —)■ eq) imposes the constraint 




k.m 


100 GeV\ 

miM ) 


+ 




k.m 


100 GeV 

TTliM 

•"m 


< 7.9 X 10 
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while the result from the Brookhaven E821 experiment on Aa^ (Eq. (41)) imposes 

100 GeV' 


k,m 


TYliM 

•"m 


< 8.6 X 10 


-4 


In Figs. (|4|)-(|^ we plot the contour of LogiQi?(/i —eq) (upper panel) and 
Log^gAa^ (bottom panel) in the (fi'osoris, ^mirror) plane for both normal (left panel) 
and inverted (right panel) neutrino mass hierarchies for scenarios 1 (red curves) and 
2 (blue curves) with the six cases of couplings aforementioned: (1) Qqs 7^ 0, gis = 0 
(Fig. 3 , ( 2 ) gis = 10 " 2 x 5 (os (Fig.|^, (3) gis = 10"^ x 5(05 (Fig.j^, (4) gis = 0.5 x 5(95 
(Fig. 0, (5) gos = gis (Fig. |§, and (6) g^s = 0, gis 7 ^ 0 (Fig. g, respectively. 

At the upper panel of each of these hgures, the (light) gray area is excluded by the 
current limit of Log]^Qi?(/i —?• eq) = —12.24 from MEG experiment [20] for scenario 
(1) 2 respectively. The projected sensitivity of Log]^Qi?(/r —>■ eq) = —13.40 [21] is also 
shown for each scenario in the two plots in the upper panel for comparison. 
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LogjQB(//—»ey) Scenario 1/2 (Normal) 



-5 -4 -3 -2 -1 0 


^ogigos) (with §15=0) 


LogjQB(/i—> 67 ) Scenario 1/2 (Inverted) 



-5 -4 -3 -2 -1 0 

Log( 5 os) (with g, 5=0) 


Log^pAa^ Scenario 2/1 (Normal) LogjQAa^ Scenario 2/1 (Inverted) 




FIG. 4. Contour plots of LogiQB{fi —> e'y) (top panel) and Log^gAa^ (bottom panel) 
on the (^ 05 )-^mirror) plane for normal (left panel) and inverted (right panel) hierarchy in 
scenarios 1 (red curves) and 2 (blue curves) with gos = gQg and gis = 9[s = 0- For details 
of other input parameters, one can refer to the text in section VI. 

At the bottom panel of each of these hgures, the red (blue) area is dehned by the 
Log^gAa^ = —8.54 |22] from the E821 experiment of the Brookhaven National Lab 
(BNL) for the discrepancy between the SM model prediction and the measurement 
for the muon anomalous magnetic dipole moment for scenario 1 (2), respectively. 

From all the plots in these hgures, we observe the following general features. 
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LogjQB(//—»ey) Scenario 1/2 (Normal) 



-5 -4 -3 -2 -1 0 


Log(g05)(withg,5=10 ^-gos) 


LogjQB(/i—>ey) Scenario 1/2 (Inverted) 



-5 -4 -3 -2 -1 0 

Log(g05) (with gi s= 10"^-gos) 


LogjQAa^ Scenario 2/1 (Normal) 


LogjQAa^ Scenario 2/1 (Inverted) 



FIG. 5. Same as Fig. S with gos = 5os and gis = g[s = 10 ^ • gos instead. 

• In the same mass range of the mirror leptons the LFV process /i —)■ ey is more 
sensitive to the couplings by almost two order of magnitudes as compared with 
the anomalous magnetic dipole moment of the muon. This is partly due to 
the fact that the i?(/r —>■ ey) is quartic in the couplings, while in Aa^ they are 
quadratic. 

• As one turns on the triplet coupling gis from 0 to gis = dos (Fig- (|^ to 
Fig. ([^), the contours for LogiQB(n —)■ ey) (upper panels) are shifting toward 
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LogjQB(//—»ey) Scenario 1/2 (Normal) 



-5 -4 -3 -2 -1 0 

LogC^os) (with gis=10"‘-go5) 


LogjQB(/i—>ey) Scenario 1/2 (Inverted) 



-5 -4 -3 -2 -1 0 

Log(g 05 ) (with gi s= 10"‘ -gos) 


LogjQAa^ Scenario 2/1 (Normal) LogjQAa^ Scenario 2/1 (Inverted) 



FIG. 6. Same as Fig. Q with gQs = g^g and gig = g[g = 10 ^ • gog instead. 

to the left, indicating the role of the triplet singlets become more relevant 
and thus the constraints on parameter space become more stringent from the 
current MEG limit. However in the last case of Fig. (|^ when the singlet 
coupling gog is set to zero such that only the triplet singlets are contributing 
in the loop diagram, the contours of LogiQB{fi —)■ ey) are slightly shifting 
back toward to the right. Similar behaviors can be found for the contours of 
Logj^o-Ao/i, but the effects are tiny and not easily seen on the log scale, except 
for the last three cases of Figs. 0-0 (lower panels). 
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LogjQB(//—»ey) Scenario 1/2 (Normal) 



-5 -4 -3 -2 -1 0 

Log(§os) (with gi i=0.5-go5) 


LogjQB(/i—> 67 ) Scenario 1/2 (Inverted) 



-5 -4 -3 -2 -1 0 

Log(gos) (with gi 5=0.5-gos) 


LogjpAa^ Scenario 2/1 (Normal) LogjQAa^ Scenario 2/1 (Inverted) 



FIG. 7. Same as Fig. S with gos = g'os and gis = g[s = 0-5 • gos instead. 

Regarding the sensitivity on the two scenarios, we can obtain the following state¬ 
ment by comparing the red and blue contours corresponding to the scenarios 1 and 
and 2 in each of these hgures. 

• The sensitivity of the couplings in the i?(/r —> ey) has been weakened by one 
to two order of magnitudes for scenario 2 as compared to scenario 1. This is 
due to the fact that in scenario 2, the three unknown unitary mixing matrices 
are now departure from which allows the couplings take on larger val- 
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LogjQB(//—»ey) Scenario 1/2 (Normal) 



-5 -4 -3 -2 -1 0 


Log(§os) (withg,5=go5) 


LogjQB(/i—>ey) Scenario 1/2 (Inverted) 

800 
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> 

o 

^ 600 
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i 

? 500 

I 

I 400 

E. 

S 300 

i 

200 


100 

-5 -4 -3 -2 -1 0 

Log(gos) (withg,5=gos) 



LogjpAa^ Scenario 2/1 (Normal) LogjQAa^ Scenario 2/1 (Inverted) 



FIG. 8. Same as Fig. Q with gos = Qqs = 9iS = 9is instead. 


ues since the amplitudes involve products of the couplings and the elements of 
mixing matrices. However this sensitivity is not present for the muon anoma¬ 
lous magnetic dipole moment as the distance between the two red and blue 
contours for the two scenarios in the lower panels of all these plots are well 
within a small range of the coupling g^s (or gis in Fig. ([^). For example, at 
Mmi rrnr = 100 GeV, the allowed value of gos varies from 10“^'® to 10“^'® (10“^'® 
to 10“^'^) as seen from the upper (lower) panels of Figs.Q-([^. 
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LogjQB(/^—>ey) Scenario 1/2 (Normal) 



-5 -4 -3 -2 -1 0 


LogUis) (with go5=0) 


LogjQB(/i—>ey) Scenario 1/2 (Inverted) 
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100 

-5 -4 -3 -2 -1 0 

Log(gis) (with goi=0) 



LogjpAa^ Scenario 2/1 (Normal) LogjQAa^ Scenario 2/1 (Inverted) 



FIG. 9. Same as Fig. Q with gos = = 0 and gis = g'lg instead. 

Regarding the sensitivity on the neutrino mass hierarchies, one can obtain the 
following statements by comparing the left and right panels in each of these hgures. 

• As one slowly turns on the triplet coupling = 0 (Fig. (|^) to 5^15 = 
10“^ X gos (Fig. (|^), the red contours of Log]^oR(/r —>■ ey) of scenario 1 in the 
left and right panels in all these plots remain the same, while the blue contours 
of scenario 2 in the right panels move toward to the left. This indicates that 
noticeable differences in the contours of Log^gR {fi —)■ ey) between the normal 
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and inverted neutrino mass hierarchies can be seen in these cases. In general 
the couplings are about an order of magnitude more sensitive in the inverted 
mass hierarchy than the normal one for scenario 2. However, for gis > 0.5 x 
these differences diminish. 

• There are no discernible differences between the two mass hierarchies for the 
muon anomalous magnetic dipole moment in both scenarios for all 6 cases of 
couplings. 

VII. IMPLICATIONS 

The constraints on the Yukawa couplings coming from /r —>■ ey has several impli¬ 
cations among which two are particularly relevant. 

• The size allowed for the Yukawa couplings by present limits on H(/i —)■ ey) 
has an important implication on the decay lengths of the mirror leptons. It 
is beyond the scope of this paper to discuss this in detail here but a few 
remarks are in order. In the search for mirror leptons, one would like to look 
for characteristic signatures which can be distinguished from SM background. 
One of such signatures could be events with displaced vertices, in particular 
events with decay lengths which are macroscopic (/ > 1mm). How this type 
of events can be correlated to /r —>■ ey is a topic which was already mentioned 
in [3]. With the present update which includes a more detailed analysis taking 
into account mixings in the lepton sector, one can have a better idea of the 
correlation between the feasibility to observe /x —>• ey and the detection of 
mirror leptons. 

A mirror lepton can decay directly into SM leptons with an accompanying 
Higgs singlet. For example, one can have Ihj + 4>ks where i,j = e, p, r 
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and k = 0,1, 2, 3. The decay length will depend on the magnitude of the 
Yukawa couplings as well as on the various mixing parameters contained in 


Eq. (19). We just take one example here for the sake of discussion. The 
interaction Lagrangian for —)■ lij + 4>ks can be expressed as { 61^12 + 

P‘L-M .22 + TlM 32 )where (for scenario 2 with the normal hierarchy) 


Mu = (5.834 X 10"® - O.OOOO25i)^os0os + (42) 

(^ 15 ( 0.324 + 0.159i) + ^25(0.407 - + 

(^is(0.154 + 0.200i) + ^2 s( 0.192 + O.238f))02S + 

(<715(0.074 - 0.325*) + <725(0.201 - 0.102*))</.35 

■M 22 = O.9999335<o50os + 

(^i5(-0.262 + 0.332*) + ^25(-0.262 - 0.332i))</)i5 + 
(^15(0.146 - 0.1930 + ^25(0.146 + 0.1930)025 + 

(^15(0.067 - 0.2550 + ^25(0.067 + 0.2550)^35 
J^32 — (0.00006 T 0.0000205^05005 T 

(5i5(-0.054 - 0.2760 + 525(-0.145 + 0.2570 ) 0i5 + 
(5i5(-0.163 - 0.0430 + ^25(0.269 + 0.405*))025 + 

(<715(0.166 - 0.503*) + <725(-0.157 - 0.077*))035 


Depending on the particular search (e, /* or r), a displaced vertex might occur. 
For instance, if one focuses on r, and if < 7^5 <C 505 , the constraint on <705 < 10“^ 
(see the above hgures) implies that 7 *^ —)• + 0 a :5 would have a macroscopic 
decay length. There are many such cases but it is beyond the scope of this paper 
to discuss this issue at length. We merely point out the relationship between 
the constraints coming from /* —)■ ey and the implication on the search for 
mirror leptons. 
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• The other implication concerns the VEV of the singlet Higgs helds. Since 
the seesaw mechanism implies the masses of the light neutrinos are given by ~ 
mjy/M and with M ~ 0{Aew), if was stated in [T] that mu ~ 0(100 keV) and 
that the singlet VEV ~ 0(100 keV) if gs ~ 0(1). However, constraints from 
/i —>• ey imply qqs < 10“^ which now brings the singlet VEV up to 0(100 MeV). 
In fact it can even be of the order 0(1 GeV). From this observation, it is safe 
to say that there does not appear to be much of a hierarchy problem in the 
EW-scale vr model. 


VIII. CONCLUSIONS 

In this work, we present an update on a previous analysis [3] for the process 
p —?■ ey performed in the original EW-scale vr model [1] to an extended model [ 6 ]. 
Mixings effects of neutrinos and charged leptons constructed with a H 4 symmetry as 
recently studied in [7] are also taken into account. In this context, the rare process 
p —)■ ey is link to interesting new physics beyond the SM in the lepton sector, like 
neutrino and charged lepton mass mixings, neutrino mass hierarchies, mirror leptons 
as well as singlet and triplet scalars of H 4 , etc. The related muon anomalous magnetic 
dipole moment is also studied in detail for the model. 

To summarize, we hnd that 

• One can deduce more stringent constraints on the parameter space of the EW- 
scale fR model by using the LEV process p —ey than the muon anomalous 
magnetic dipole moment. 

• The branching ratio i?(/i —)■ ey) shows some sensitivity to the neutrino mass 
hierarchies in scenario 2 but not scenario 1 , depending on the ^4 triplet cou¬ 
pling constants. However we are not advocating the use of the process p —>■ ey 
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to settle the issue of neutrino mass hierarchies. After all, this is a rare process. 


• More stringent constraints can be deduced in scenario 1 than scenario 2 using 

B(/i ^ 67). 

• Future data from MEG experiment with the projected sensitivity will impose 
further constraints on the parameter space of the model. 

• The muon anomalous magnetic dipole moment is sensitive neither to the neu¬ 
trino mass hierarchies nor the scenarios for all 6 cases of the couplings studied 
here for the model. 

Searching for new physics via rare processes is complementary to direct production 
of new particles at colliders. For fi —)■ ey, the relevant new particles in the model are 
the mirror leptons and scalar singlets running inside the loop diagram. As shown 
in our analysis, the Yukawa couplings of the Higgs singlets to the leptons in the 
EW-scale model are constrained to be small in order to be consistent with the 
current experimental limit on i?(p —)■ ey). Thus searching for mirror particles of 
this model at the LHC would be quite interesting since, due to small couplings, they 
might decay outside the beam pipe and inside the silicon vertex detectors. The A 4 
singlet and triplet scalars are likely to escape detection as missing energies. 

As an outlook, one would like to generalize this work to fi — e conversion. This 
work is now in progress and will be reported elsewhere [23]. 


APPENDIX 


For the general case of retaining the external fermion masses rriij, the integrals 


•PkS 


and J' 


•t-kS 


m 


Eqs. (29 )-(30) have to be replaced by X j 


JIL 1 

<PkS 

\ 

\ 
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and J ( ^ ) respectively, where 

\ "SM / 

\ ‘'Tn ‘'Tn '■m / 

f\ r^y -— . 

Jo Jo a; + 1 /+ (1 - a; - |/)(r - xrj - yrd - *0+ 

Jo Jo ^^x + y + {l-x-y){r-xrj-yri)-iO+ ' 
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